Photoemission Spectroscopy from Inhomogeneous Models of Cuprates 

J. Eroles^'2, G. Ortiz\ A. V. Balatsky^ and A. R. Bishop^ 
^Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. 
^Centro Atomico Bariloche and Instituto Balseiro, S. C. de Bariloche, Argentina. 
(Received February 1, 2008) 

We investigate the electronic dynamics in the underdoped cuprates focusing on the effects of 
one- dimensional charge stripes. We address recent experimental Angular- Resolved Photoemission 
Spectra results on (Lai.28Ndo.6Sro.i2)Cu04. We find that various inhomogeneous models can ac- 
count for the distribution of quasiparticle weights close to momentum k = (7r,0) and symmetry 
related points. The observed flat dispersion region around the same k point can only be addressed 
by certain classes of those inhomogeneous models which locally break spin symmetry. Homogeneous 
models including hopping elements up to second neighbors cannot reproduce the experimental quasi- 
particle weight, since most of it is centered around k = ^). 



Recent interest in the inhomogeneous states of high- 
Tc superconductors has gained further impetus from new 
experiments that support charged stripe formation. One 
set of data is available on the Nd-doped LSCO [Q and, 
more recently, a second set concerns oxygen doped LSCO 
[||. Angular- Resolved Photoemission Spectra (ARPES) 
obtained on the commensurate-doped Nd LSCO [|j re- 
vealed that the electronic spectra, believed to contain 
static stripes, are consistent with the electronic states 
having quasi-one-dimensional (ID) character. 

The question that naturally arises from these measure- 
ments is: what are the most appropriate models that can 
capture both the formation of the stripes and adequately 
describe the spectroscopic features of the electronic states 
in these compounds? The issue of the physics that drives 
stripe formation is a matter of extensive current debate 
[H , with a number of models that can describe the static 
and dynamical stripe formation. 

Faced with the new experimental data, we choose here 
a different approach. We assume from the beginning 
that the Hamiltonian for our model is inhomogeneous 
to mimic the effects of stripes. As a consequence the 
electronic states will be quite different from the conven- 
tional homogeneous models, such as a tt'-J model. We 
find that suitable explicitly inhomogeneous models allow 
a good fit to the ARPES data, and capture the main 
qualitative features that are consistent with a quasi- ID 
nature of electronic states on the stripes. 

The inhomogeneous models we will consider have al- 
ready been introduced in Ref. ^ , where the basic micro- 
scopic scenario starts from a homogeneous t-J Hamilto- 
nian as a reference background model: 

Ht-J = -t + J ^(Sr • Sr - ^nrnr). (1) 

(r,f), cr (r,r) 

We have also considered the tt'-J model which includes 
hopping elements up to second neighbors. To mimic the 
stripe segments, we add specifically inhomogeneous mag- 
netic interactions. These inhomogeneous terms locally 



break translational invariance and spin-rotational SU (2) 
symmetry: 

H,^^=Y,6J,S^^SI^^-^[S^S^^S-S;). (2) 

Here 6J±_ 7^ 6Jz, represents the magnetic perturbation 
of a static, local, Ising anisotropy, locally lowering spin 
symmetry (termed a t-JJz model Only a few links 
(where the stripes are located) have this lowered spin 
symmetry. The Ising anisotropy is also sufficient to pro- 
duce a spin gap ||] and pair binding of holes in this class 
of models is substantial [||. An alternative scenario for 
stripe segments can be realized through inhomogeneous 
terms of the form 

^inh = -^6,(l-n,) , (3) 

i 

where > represents an on-site orbital energy which 
is non- vanishing only at the stripe sites i {tt'-Jci model). 
ARPES signal of the noninteracting electron model with 
external potential, consistent with the stripe pattern, has 
been previously addressed in Ref. [^. The spin rotation 
symmetry is not broken in this case. This kind of model 
resembles, although is different, to that studied in Ref. 

, where it was shown to reproduce experimental neu- 
tron magnetic scattering, optical conductivity and the 
ARPES measured broad spectrum near k = (f , f ) along 
the direction (0,0) - (tt, tt) 0. 

In the rest of this paper the energy scale will be deter- 
mined by t. As we will vary the magnitude of the spin 
anisotropy we choose to fix J/t = 1. This is a little 
higher than the commonly assumed value for homoge- 
neous cuprate systems of J/t = 0.4, but close enough to 
be physical. It is not our purpose here to define an energy 
scale which can be related to the experiments, but rather 
to demonstrate that the main aspects of the experiments 
can be address by particular kinds of models. Systematic 
tuning of the parameters will be published elsewhere. 
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Recent ARPES experiments in Lai.28Ndo.6Sro.i2Cu04 
seem to indicate a ID electronic structure At low 
transferred energies and integrating the spectral weight 
within a 100 meV window {Aco =100 meV), the quasi- 
particle weight looks like that depicted in Figs. and 
||. Basically, there is almost no weight at k = (f,f), 
and most of the weight is concentrated near the point 
k = (7r,0). Moreover, one can see from the quasiparti- 
cle band structure that there is almost no dispersion in 
ky. Besides, there is a flat band for large kx that can be 
related to a velocity scale By integrating the same 
data over a window of 500 meV {Aco =500 meV), the 
experimental spectrum is like that in Fig. ^. There is no 
peak at k = (f , f ) and most of the spectral weight is 
transferred from k = (tt, 0) to k = (f ,0) and symmetri- 
cally equivalent points. Both cases {Auj = 100 meV or 
500 meV) are characteristic of a ID behavior of the elec- 
tronic degrees of freedom . These experiments were all 
performed at a commensurate hole filling corresponding 
to X = 1/8. 



where I^J^) corresponds to an eigenstate of the Hamil- 
tonian in the subspace of u holes with energy E^. The 
first poles of the corresponding Green's function define 
the quasiparticle weight through the intensity, and the 
energy associated to that quasiparticle through the posi- 
tion of the poles. As the experiments in Ref. |||] are at 
X = 1/8, we carried out the simulations with one hole 
added to a background where there was already one hole 
present {u = 1). For our 16 sites cluster, the final con- 
centration then corresponds to x = 1/8. 

Since our system is finite, we can only consider a fi- 
nite set of k values in the Brillouin zone. These are: 
k = {kx^ ky) with kx^y = 0, ±7r. To compare with ex- 
perimental results we associate a Gaussian with a small 
width (cr c^i 0.1) centered at each computed k point and 
with a height corresponding to the spectral weight for 
that k point. Finally, we sum all these Gaussians cen- 
tered at the k values of our finite lattice. The resulting 
field is color coded with blue corresponding to maxima 
and red to minima. For each graph the maximum is 
scaled to 1. 




FIG. 1. (a) Experimental ARPES from ||. (b), (c) and 
(d) correspond to the unsymmetrized calculated spectrum for 
t-JJz, tt'-J and tt'-Jti models, respectively. Note that no 
integral of spectral weight has been taken yet. 

Our goal here is to study different homogeneous and 
inhomogeneous models and compare to the experimen- 
tal dynamic properties measured in Ref. We used 
the Lanczos numerical exact diagonalization method on 
a 4x4 cluster (see inset Fig. ^). When stripes are present, 
the Hamiltonian along one leg is modified by adding the 
inhomogeneous terms (like in Eqs. ^ or ||). 

The spectral function we need to compute is given by 

^(k,c.) = ^ |(*r'|cLl*o)l' 5{u - (iS^i - K)) , 

n 

(4) 
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FIG. 2. Comparison between the inhomogeneous t-JJz 
model (t/J = 1 and 5J± = —0.9) integrated up to Ao; = 0.13t 
(left) and ARPES experimental results (right) from in- 
tegrated over Aix) — lOOmeV. (a) Four quadrants view (b). 
Zoom of the zone around (tt, 0). 

In Fig. we show (a) the experimental spectrum from 
Ref. [|| and the corresponding unsymmetrized model pre- 
dictions for (b) t-JJz, (c) tt'-J and (d) tt'-Jci models. 
For the calculations we have taken the first quasiparti- 
cle pole for each k point. (Note that for some k points 
the weight is very small). Clearly predictions from the 



2 



inhomogeneous models ((b) and (c)) are qualitatively 
different from the homogeneous tt'-J model. In both 
inhomogneous cases the weight is concentrated around 
k = (tt, 0) as in the experimental ARPES. But if stripes 
are present, the experiment averages domains with dif- 
ferent stripe orientations. Therefore, in order to compare 
with the experiments, we symmetrize the numerical re- 
sults by taking the average with the rotated spectrum 
and integrating the results in a frequency window Acj 
for each k point, as was done in Ref. P]. 
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FIG. 3. Spectra for different parameters and integra- 
tion windows for the t-JJz model. As the strength of the 
stripe is decreased {3J± 0, making the model more homo- 
geneous) noticeable weight is transfered to k = (f^f) and 
to the boundary of the "ghost" Anti-Ferromagnetic Brillouin 



FIG. 4. Comparison between the inhomogeneous t-JJz 
model (t/J = 1 and 5Ja_ = -0.9) (left) and ARPES experi- 
mental results (right) [|3| integrated in a 500 meV range. Four 
quadrants view. (a)Au; = 2.5t; (b)Au; = 3.5t. 

In Fig. ^ we compare the calculated spectra in the in- 
homogeneous t-JJz model (Eq. ||), where the stripes are 
represented as local Ising anisotropy, symmetrized and 
integrated over a range up to Ao; =100 meV. The Hamil- 
tonian corresponds to Eqs. and ^ with J Jx = —0.9. Fig- 
ure ^ corresponds to (a) the full 4 quadrants view, and 
(b) to a partial view around k = (tt, 0). The agreement is 
very good. In the experiment significant spectral weight 
was found around k = (0,0). Our data, on the other 
hand, do not show any significant weight at k = (0,0). 
It is well known that the reduction to one band effective 
models is not very good at this k point [^. Thus the 
lack of weight at k = (0, 0) is a consequence of the reduc- 
tion to the t-J model and not important for our present 
considerations. 

In Fig. H we display the quasiparticle weights for var- 
ious values of the stripe anisotropy. Note the transfer 
of weight from k = (tt, 0) to k = (f , f ) as the system 
becomes more homogeneous. 

Quasiparticle features involving many sites (and there- 
fore with higher energies) will not be properly captured 
in our small cluster. Nevertheless, it still seems able to 
describe the main trends as the window of integration in 
the experiment is raised to Ao; = 500 meV. As seen in 
Fig. |^(a) the weight concentrates around k = (|-,0) and 
equivalent points, even though the weight at k = (f , f ) 
is still significant. As the window of integration is raised 
to Acj = 3.5t, most of the weight is concentrated between 
k = (0,0) and k = (f ,0). Probably because of the size 
of our cluster, none of the models studied can exactly 
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In Fig. I we compare the distribution of quasiparticle 
weights for homogeneous tt'-J and inhomogeneous t-JJz 
models after symmetrization and for different integration 
windows Acj. We study the tt'-J model since it is well 
known that, in our cluster, the bare t-J has a hidden sym- 
metry for 1 hole that makes k = (tt, 0) and k = (f , f ) 
degenerate. In particular for a chain, the dispersion re- 
lation for a t' term is Ck = — cos(2/c), and therefore it 
lowers the quasiparticle energy at both k = and, more 
importantly, at /c = tt. If for some reason the system has 
a ID (or quasi ID) behavior, this f term will transfer 
quasiparticle weight from /c = |- to /c = tt, as observed in 
the experiments. From Fig. || it can be seen that in the 
homogeneous model the experimental weights cannot be 
reproduced by t' /t = —0.3 (the commonly assumed value 
for f [^) for any Aco. It is also interesting to note that in 
the inhomogeneous tt'-Jei model, after symmetrization, 
the prominent peak at k = (tt, 0) and equivalent points 
(as in Fig. ^(b) for Acj = O.Olt and Acj = O.lt) is com- 
parable to k = (^, for most of the Ao; window sizes. 
Therefore the agreement with experiments is poor. For 
large Acj all the models are very much alike. This is to be 
expected since if the stripe perturbation is small enough, 
only the low energy states change significantly, and there- 
fore all inhomogeneous and homogeneous models give the 
same qualitatively spectra for sufficiently large Acj. 
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FIG. 5. Comparison of the spectral density between (a) 
homogeneous tt'-J {t' = — 0.3t) and (b) inhomogeneous t-JJz 
{t/ J —1 and 5J± — —0.9) models. Blue corresponds to max- 
imum intensity and red to minimum. Each plot corresponds 
to the indicated range of integration Auo (in units of t). 

As mentioned earlier, another distinctive feature in 
the ARPES experiments is the ffat dispersion around 
k = (7r,0). This ffat region may be associated with a 
velocity which in turn defines the linear relation between 
the separation between stripes and the superconducting 
transition temperature [^. Therefore it is important to 
see if the fiat region is present in our models. In Fig. ^ we 
show the position of the first pole for each model, for the 
k points along the (0, 0) - (f, 0) - (tt, 0) line. The t-JJ^ 
and tt'-J models are fiat enough to be consistent with the 
experiment, but as we have seen above, the latter cannot 
account for the weight distribution (see Fig. ||). On the 
other hand the model labeled as tt'-Jci does not show 
the fiat dispersion. This is a further indication that even 
though inhomogeneities are necessary to account for the 
ARPES experiments, not all kinds of inhomogeneities are 
appropiate. 
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FIG. 6. Dispersion relation for the different models along 
the (0,0) — (f ,0) — (tt, 0) line. The experimentally observed 
flat band close to k = (tt, 0) |^ is well reproduced by the 
t-JJz and tt'-J models. 

In conclusion, we have calculated the ARPES response 
in several inhomogeneous and homogeneous models. We 
used exact diagonalization in a cluster with C/^y symme- 
try and hole doping appropiate to x = | . We found that 
the experimentally measured quasiparticle peak near k = 
(tt, 0) and symmetry equivalent points is well reproduced 
only by the inhomogeneous models. We also find that 
for all models, some weight peaked at k = (^, ^). It is 
of comparable intensity (although usually smaller) to the 
peak at k = (7r,0), and is present for some integration 
windows. Of those models we studied, only those that 
locally break SU(2) symmetry can also explain the ex- 
perimentally observed flat dispersion around this same 
k point. We note that this same model has a spin-gap. 



substantial binding of holes and exhibits an incommen- 
suration in the spin structure factor It thus ap- 

pears to capture many important experimental features 
in the electronic, magnetic and potentially superconduct- 
ing channels. 

Preliminary results of this work were presented on 
the M2S-HTSC-IV conference in Houston Subse- 
quent to completion of this wor k we learn about a r ecent 
preprint by D. Orgad e t al. (|cond-mat/0005457| ) and 
by M.G. Zacher et al ( |cond-rnat/0005473D , where the 
relationship between stripe patterns and ARPES data 
is discussed. We thank A. A. Aligia, S. Kivelson, Z-X. 
Shen and X.J. Zhou for useful discussions. Work at Los 
Alamos is sponsored by the US DOE under contract W- 
7405-ENG-36. 
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